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ABSTRACT 

For  fiber  reinforced  cement  based  composites,  the  principal  beneficial 
effects  of  fibers  (metallic,  mineral  or  organic)  accrue  after  the  matrix 
has  cracked.  For  loads  beyond  which  the  matrix  has  initially  cracked,  the 
further  crack  extension  and  opening  is  resisted  by  bridging  of  fibers  across 
the  cracks.  The  resistance  provided  by  the  fibers  will  depend  principally 
on  the  debonding  and  the  pull-out  resistance  of  fibers. 

A  theoretical  model  based  on  the  concepts  of  nonlinear  fracture  mechanics 
to  predict  the  resistance  provided  by  the  fibers  against  the  fracture  of 
matrix  is  presented  in  this  report.  The  theoretically  predicted  response 
is  compared  favorably  with  the  experimental  data  on  notched  beams  and  double 
cantilever  beam  specimens  of  steel  fiber  reinforced  concrete.  The  proposed 
theoretical  model  provides  a  method  to  calculate  fracture  resistance  for  a 
crack  extension  in  a  specimen  of  any  geometry. 

One  of  the  key  parameters  required  for  the  model  is  the  relationship 
between  the  uniaxial  post-cracking  stress  and  the  corresponding  displacement. 
This  relationship  will  depend  on  the  bond-slip  function  of  fibers.  A  method 
to  estimate  this  relationship  is  presented*. 
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INTRODUCTION  AND  SUMMARY  OF  CONCLUSIONS 


Research  conducted  during  the  last  twenty  years  has  shown  that  the 
addition  of  fibers  significantly  Improves  penetration,  scabbing  and  frag¬ 
mentation  resistance  of  concrete.  The  possible  applications  of  fiber 
reinforced  concrete  (FRC)  include  explosion  and  shock  resistant  protective 
structures.  Even  though  the  enhanced  "cracking  resistance"  is  the  most 
important  attribute  of  FRC,  there  are  no  rational  methods  of  measuring  or 
predicting  this  important  materialproperty . 

For  fiber  reinforced  cement  based  composites,  the  pricipal  benefi¬ 
cial  effects  of  fibers  (metallic,  mineral  or  organic)  accrue  after  the 
cracking  of  matrix  has  occured.  For  loads  beyond  which  the  matrix  has 
initially  cracked,  the  further  crack  extension  and  opening  is  resisted 
by  bridging  of  fibers  across  the  crack.  To  incorporate  the  effects  of 
fibers  bridging,  many  investigators  have  used  the  classical  linear  elastic 
fracture  mechanics  concepts.  These  past  attempts  however  assume  one  or 
more  of  the  following:  (1)  the  closing  pressure  due  to  the  fiber-bridging 
is  constant,  (2)  the  extent  of  the  fiber  bridging  zone  is  small  compared 
to  the  length  of  the  traction  free  crack  length,  and  (3)  ignore  the  energy 
absorbed  in  the  fiber-bridging  zone.  In  this  report  the  results  of  a 
theoretical  model  which  does  not  make  the  above  assumptions  are  presented. 

An  existing  crack  in  a  cement  based  matrix  is  replaced  by  an  effective 
crack  length  which  consists  of  three  zone:  (1)  a  traction  free  zone, 

(2)  fiber  bridging  zone  and  (3)  the  matrix  process  zone  (E^) 

resulting  from  the  aggregate  interlock  and  microcracking  (Fig.  5).  It  is 
assumed  that  the  stresses  and  deformations  due  to  the  applied  loads  and 


Che  closing  pressures  can  be  calculated  using  theory  of  elasticity.  The 
fiber  bridging  closing  pressure  depends  on  the  crack  opening  displacement 
which  in  turn  depends  on  the  geometry  of  the  specimen,  external  loading 
and  on  the  closing  pressure  itself.  An  iterative  procedure  was  developed 
to  account  for  this  coupling  between  the  closing  pressure  and  the  crack 
opening  displacement . 

Experiments  with  fiber  reinforced  concrete  specimens  subjected  to  uni¬ 
axial  tension  in  a  relatively  stiff  testing  system  indicate  that  the  post¬ 
peak  displacement  are  essentially  a  result  of  opening  of  a  single  crack. 

Thus  it  was  assumed  that  the  relationship  between  fiber  bridging  closing 
pressure  and  the  crack  opening  displacement  (COD)  are  equivalent  to  the 
uniaxial,  tensile,  post-peak  stress-displacement  function.  A  method  to 
estimate  this  relationship  is  proposed  (Fig.  8)  for  steel  fiber  reinforced 
concrete.  Experiments  to  verify  this  concept  are  currently  underway  at 
Northwestern  University.  Based  on  these  experiments  it  is  expected  to 
derive  closing  pressure  vs.  crack  opening  displacement  relationship  from  a 
single  fiber  pull-out  vs.  slip  function. 

Based  on  the  theoretical  model,  load  vs.  crack  mouth  displacement  curves 
for  notched,  steel-fiber  reinforced  concrete  beams  were  derived  and  compared 
with  experimental  data  (Fig.  26).  The  crack  opening  displacements  for  the 
theoretical  model  were  calculated  using  some  simple  approximations.  A  more 
rigorous  method  of  calculating  COD  by  modelling  crack  as  distributed  dislo¬ 
cations  is  currently  underway. 

Fracture  resistance  vs.  crack  extension  relationships  (R-curves)  were 
calculated  and  compared  with  the  experimental  data  on  the  double  cantilever 
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beam  specimens  of  steel  FRC  (Fig.  20).  The  theoretical  curves  were  obtained 
by  assuming  that  the  displacement  resulting  from  the  fiber  bridging  forces 
are  irreversible.  A  good  correlation  was  obtained  with  the  experimentally 
observed  R-curves  as  well  as  load  vs.  load-line  displacement  curves  (Fig. 22). 
Note  that  if  the  asymptotic  value  of  R-curves  is  considered  a  material 
parameter,  then  that  value  can  be  a  useful  quantity  in  identifying  the  bene¬ 
fits  of  fiber  addition.  For  example  the  predicted  steady  state  value  for 
specimens  reinforced  with  37.5  mm  steel  fibers  is  approximately  40  times 
that  for  plain  unreinforced  matrix  (Fig.  25).  This  relative  improvement  in 
fracture  energy  is  comparable  to  that  obtained  by  using  an  empirical  method 
called  toughness  index  suggested  by  American  Concrete  Institute. 
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PREVIOUS  THEORETICAL  MODELS 


Many  Investigators  have  attempted  to  predict  the  fracture  behavior  of 
fiber  reinforced  composites  based  on  the  fiber  bridging  concept  £2,  3,  4, 

12,  13]. 

Lenian  and  Bunsell  £2]  proposed  that  the  stress  intensity  factor  (Kjj) 
required  for  crack  propagation  was  the  algebraic  sum  of  the  stress  intensity 
factors  resulting  from  the  external  load  applied  to  the  cracked  specimen  and 
without  considering  the  presence  of  fibers  (Kq)  and  the  stress  intensity 
factor  due  to  fibers  bridging  across  the  crack  (1^) ,  i.e., 

Kr  -  Kq  -  Kr  -  Kc  (1) 

where  Kc  *  experimentally  determined  critical  stress  intensity  factor  for 
the  unreinforced  matrix. 

They  calculated  the  effect  of  fibers  bridging  (Kr)  in  Asbestos  Cement 
by  employing  Paris  and  Sih's  analytical  expression  £l4]  which  was  derived 
for  an  Internal  crack  in  an  infinite  sheet  under  a  region  of  uniform  closing 
pressure  9  applied  at  the  crack  tip  (Fig.  1-a) .  This  relationship  is  given 


as: 


^  *  'Ipr^in'1  f  -  sin’1  |  -  (1  -  +  (1  -  j 


(2) 


where  a,  b,  and  c  are  crack  length,  and  the  boundaries  of  the  applied  uni¬ 
form  pressure  9  ,  as  shown  in  Fig.  l-a. 

They  assumed  that  the  closing  stresses  produced  by  the  asbestos  fibers 
were  the  same  for  microcracks  as  for  the  main  crack,  and  that  a  zone  of 
microcracking  of  length  can  be  represented  by  an  imaginary  increase  of 
crack  length,  xjj^,  where  x  is  a  factor  less  than  unity,  to  which  the  closing 
stress  was  uniformly  applied  (Fig.  1-b). 
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The  stress  intensity  factor  (K0)  was  then  calculated  as: 


**  "  K0  -y=U  +  X*P) 


4  /  Tr 


2  -  sin 


-1  - v  +  [l  + - ^ - -1  h) 

(a  +  x*p)  |_  (a  +  xJLp)2J  / 


where  x  was  the  fraction  of  the  microcracked  zone  J,p.  Note  that  a  is  the 
uniform  bridging  pressure  of  asbestos  fibers  which  was  assumed  equal  to  the 
maximum  value  of  the  post-cracking  stress  (*max)  and  was  calculated  from 

[15]: 


A1  Vm  +  4  \  "e  T  Vf  -a 


(A-a) 


Where  is  a  factor  related  to  the  extent  of  matrix  cracking  and  which  is 

equal  to  or  less  than  unity,  V  and  o  are  the  percent  volume  fraction  of 

m  in 

the  matrix  and  the  matrix  tensile  strength.  The  second  term  in  Eq .  (A-a)  is 

primarily  the  contribution  from  fibers  being  pulled  out.  n  is  the  effi- 

l 

ciency  factor  for  the  embedment  length  of  fibers.  Pg  is  the  coefficient  for 

fiber  orientation  and  distribution,  while  r  ,  Vf,  and  -4-  are  the  bond 

d 

strength  of  fiber-matrix  interface,  fiber  volume  fraction,  and  the  aspect 
ratio,  respectively. 


The  matrix-cracking  factor  is  usually  small  and  can  normally  be  neg¬ 
lected  (equals  zero) .  The  values  of  n  ,  Va  and  r  are  generally  a  func- 

l  9 

tion  of  specimen  geometries,  casting  procedure,  type  and  size  of  fibers, 
fiber  volume  fraction,  and  fiber  distribution  (1-D,  2-D,  or  3-D). 


By  introducing  A2  *  A  n  r  r,  Eq.  (A-a)  can  be  rewritten  as: 


-  A1  Vm  %  +A2Vft 


(A-b) 


The  values  of  Aj  and  A2  are  varied  depending  on  different  testing  con¬ 
figurations,  conditions  of  matrix  cracking  and  fiber  debonding.  In  this 
study,  the  value  of  A}  was  assumed  to  be  zero  while  different  values  of  A2 
were  selected,  depending  from  different  testing  and  cracking  conditions. 
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Their  theoretical  model  did  not  provide  good  correlation  with  the  ex¬ 
periment  results.  This  may  be  explained  from  the  fact  that,  firstly,  they 
used  a  solution  which  was  derived  for  a  crack  in  an  infinite  sheet  to  rep¬ 
resent  a  specimen  with  finite  boundaries.  Secondly,  the  actual  fiber 
bridging  pressure  depends  on  the  crack  surface  displacement.  This  implies 
that  the  assumption  of  constant  fiber  bridging  pressure  is  not  an  accurate 
one.  Thirdly,  their  method  of  including  the  effect  of  the  matrix  process 
zone  by  a  factor  x  may  not  be  very  accurate  because  x  reduces  the  effective 
crack  length  which  influences  Kq  and  Ky. 

Foote,  Cotterell  and  Mai  [3]  proposed  a  model  using  the  LEFM  analyses 
for  a  semi-infinite  crack  (Fig.  2).  Thev  also  assumed  that  the  stress  in¬ 
tensity  factor  (KR)  for  propagating  a  crack  in  fiber  reinforced  concrete 
specimen  consisted  of  two  components,  i.e., 

Kr  -  K0  -  Kr  (5) 

where  K0  is  the  stress  intensity  factor  at  the  tip  of  the  crack  if  there 
were  no  fibers  bridging  and  Kr  is  the  stress  intensity  factor  due  to  the 
bridging  fibers  which  tends  to  close  the  crack  and  was  therefore  considered 
as  negative  (this  is  similar  to  Lenian  and  Bunsell's  model  [2]). 

However,  they  realized  that  the  fiber  bridging  stresses  (crack  closing 
stresses)  were  coupled  to  the  crack-surface  displacements  which  in  turn  was 
a  function  of  the  external  loading  and  the  crack  closing  stresses.  To  solve 
this  problem,  they  assumed  a  fifth  order  polynomial  for  the  crack  closing 
stress,  v(t).  By  employing  the  semi-infinite  crack  concent  and  the  ana¬ 
lytical  expression  from  Muskhelishvili's  potential  function  [l6]  which  gives 
the  displacements  of  the  crack  faces  due  to  closing  pressure  *(t),  they 
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calculated  the  displacement  and  stress  intensitv  factor  due  to  fiber 


bridging  as: 


(6) 

(7) 

(8) 


where  ur  and  Kr  are  the  displacement  and  the  stress  intensity  factor  due  to 
the  bridging  fibers. 

The  displacement  for  the  crack  without  fiber  bridging  (uQ)  was  deter¬ 
mined  from  theory  of  elasticity  as: 


where  Poissons'  ratio  v  was  neglected  in  Eq.  (9),  (  v  *  o) . 

With  the  calculated  total  displacement,  u  •  uQ  +  ur,  an  Iterative 
method  was  adopted  to  determine  the  curve  for  asbestos-cement  composite, 
using  Che  initiation  criterion  as  equalled  Kc  (a  critical  initiation 
value  for  matrix)  : 

kR  “  Ko  +  ^  “  Kc  (10) 

The  comparison  of  the  predicted  results  with  the  experimental  values 
was  poor.  They  attributed  it  to  the  fact  that  the  assumed  crack  closing 
pressure  was  inaccurate  due  to  lack  of  experimental  data.  Besides,  the 
serai-infinite  crack  solution  was  clearly  not  the  actual  solution  for  a 
finite  specimen  and  crack  growth. 


Since  Che  crack  closing  stress  e(t)  on  Che  crack  surfaces  depends  on 
the  displacement,  u(r) ,  and  in  turn  the  displacement  u(r)  depends  on  the 
stress  function  <!>  (Z) ,  Eq.  (6)  cannot  be  solved  directly  without  knowing 
the  exact  ff(t)  and  therefore  an  iterative  method  is  used.  By  assuming  a 
fifth  order  polynomial  function  as  the  initial  value  of  a(t),  u(r)  was 
then  calculated  from  Eqs.  6  and  7,  and  the  crack  closing  stress  a(t)  was 
determined  and  Iterated. 

Using  two  types  of  test  specimens,  double  cantilever  beam  and  single 
edge  notched.  Bowling  and  Proves  [4]  explained  that  as  the  crack  propagated 
in  the  DCB  specimen  and  became  bridged  by  the  aligned  nickel  wires  (Fig.  3), 
che  measured  critical  stress  intensity  factor,  for  further  propagation. 
Increased . 

The  Increase  was  clearly  due  to  the  presence  of  wires  bridging  the 
crack  and  exerting  crack  closing  pressure  which  reduced  the  actual  stress 
intensity  factor  pertaining  to  the  matrix  crack  tip.  Thus  a  stress  inten¬ 
sity  factor  due  to  wires  bridging  (K^  must  be  subtracted  from  that  due  to 
the  applied  force  (Ka)  and  the  Initiation  criterion  was  set  as  the  net 
stress  Intensity  factor  (Ka  -  1^)  reached  the  critical  value  of  KjC,  i.e., 

*a  “  *w  "  KIC  (ID 

The  KjC  value  here  was  taken  as  that  required  for  propagating  an  un- 
brldged  crack  past  a  wire.  The  value  of  Kfl  for  the  case  of  DCB  was  cal¬ 
culated  from  an  expression  given  by  Wiederhorn  et  al  [l7]  as: 

K  -  - La -  (  3.457  +  2.315  -A-  )  (12) 

8  wW/2 

where  P  was  the  applied  load,  a  was  the  measured  crack  growth,  and  w,  b, 
and  d  were  specimen  dimensions  as  defined  in  Fig.  3. 
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To  calculate  the  contribution  of  the  bridging  wire  0^) ,  P  in  Eq .  (12) 
was  substituted  by  the  tensile  force  in  the  wire,  and  a  was  replaced  by  the 
distance  from  the  wire  to  the  crack  tip. 

The  main  argument  for  this  proposed  model  was  that  the  assumed  crack 
bridging  force  P  was  independent  of  crack  opening  disolacement  which  is  un¬ 
likely  to  be  the  case  of  fiber  reinforced  concrete  [3,  12,  13,  19,  20,  21 J. 
The  effect  of  matrix  process  zone  in  front  of  the  crack  tip  was  also  neg¬ 
lected.  The  analysis  was  primarily  designed  for  the  aligned  fibers  and 
therefore  cannot  be  used  for  the  case  of  random  distributed  fibers.  The 


formula  given  in  Eq.  (12)  was  also  based  on  the  elastic  analysis  which  is 


unlikely  to  be  applicable,  without  any  modification,  to  nonlinear  materials 
such  as  concrete  and  fiber  reinforced  composites. 

Petersson  and  Hillerborg  [l2,  13]  introduced  a  fictitious  crack  model 
using  finite  element  analysis  to  determine  fracture  toughness  of  concrete 
and  fiber  reinforced  composites.  Their  model  was  based  on  the  Dugdale- 
Barenblatt  concept  [10,  18]  which  primarily  depends  on  the  stress-  displace¬ 
ment  relationship  in  the  nonlinear  zone  (fiber  bridging  and  matrix  process 
zone) .  Their  stress-displacement  relationship  ( v  vs  ij  )  can  either  be  ob¬ 
tained  from  a  uniaxial  tensile  test  or  from  a  notched  beam  specimen.  A  set 
of  linear  stress-displacement  relationships  was  assumed  in  the  analysis. 

By  keeping  the  work  done  or  area  under  the  »  -  n  curve  constant,  the  maxi¬ 
mum  pulled-out  displacement  (  which  was  later  used  as  the  initiation 

criterlcn  (Fig.  4)  was  calculated  from: 


G 


c 


d(n  ) 


(13) 


F 
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With  the  linear  e  -  i»  relationship,  the  maximum  pulled-out  displace- 
max 

ment  <1  ^  can  be  expressed  as: 


max 

f 


(14) 


It  can  be  seen  that  the  value  of  ^max  was  calculated  as  a  function 

of  only  Gc  and  (  ffc)mav  .  However,  *»^max  should  be  primarily  a  function 
of  fiber  geometry  and  distribution.  For  example,  it  is  shown  that  the 

upoerbound  value  of  equals  whether  fibers  are  steel  or  glass 

A 

even  though  the  values  of  f/  &  may  be  different  for  these  two  types  of 

max 

fibers. 
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BASIC  CONCEPTS  OF  THE  PROPOSED  MODEL 

A  crack  jusc  prior  to  its  extension  in  Mode  I  opening  in  a  fiber  rein¬ 
forced  concrete  specimen  is  shown  in  Fig.  5.  The  length  of  the  crack  can 
be  divided  into  three  regions:  1)  a  traction-free  crack  length  which  con¬ 
sists  of  Initial  cast  notch  and  the  zone  of  crack  length  where  fibers  are 
completely  pulled  out  of  the  matrix,  2)  the  region  of  fiber  bridging  ( 
and  3)  the  matrix  process  zone  (generally  due  to  aggregate  debonding  and 
interlocking)  in  front  of  the  crack  tip.  Both  fiber  bridging  and  matrix 
process  zone  provide  resistance  to  crack  opening.  The  effect  of  the  fiber 
bridging  is  normally  much  more  significant  than  that  contributed  from  the 
matrix  process  zone,  and  as  a  result,  the  crack  closing  pressure  in  the 
matrix  process  zone  is  neglected  in  this  study. 

Region  1  (traction-free  crack)  and  region  2  (fiber  bridging)  were 
separated  at  the  point  where  the  crack  surface  disolacement  was  equal  to 
’’max  (the  maximum  crack  displacement  where  fiber  bridging  stress  is  zero, 
since  all  fibers  at  that  point  are  completely  pulled  out^  Fig.  6),  while 
the  fiber  bridging  zone  and  the  matrix  process  zone  were  divided  at  the 
point  where  crack  surface  displacement  equalled  to  (the  maximum  dis¬ 

placement  of  the  matrix  in  the  descending  branch  of  the  uniaxial  tensile 
test  where  stress  is  equal  to  zero.  Fig.  6). 

The  value  of  has  been  reported  by  Wecharatana  and  Shah  [l9,  20, 

max 

2l] .  This  value  is  about  0.8  x  10  in.  for  matrix  mix  1:2:0:0.5  (c:6:A:w). 
The  9  ^  value  used  in  this  study  was  half  the  fiber  length.  This  may  be 

justified  (as  shown  in  Fig.  6)  from  the  observation  that  fibers  are  randomly 
distributed  across  the  crack.  The  smallest  pulled-out  distance  is  equal  to 
zero  (labelled  3  in  Fig.  6)  and  the  largest  pulled-out  distance  is  half  the 


9 


fiber  length  (labelled  2  in  Fig.  6).  This  implies  that  if  the  two  crack 
surfaces  are  separated  by  a  distance  of  half  the  fiber  length,  there  will 
not  be  any  fibers  left  bridging  across  the  crack  which  subsequently  means 
the  fiber  bridging  pressure  is  zero. 

If  fibers  are  randomly  distributed  rather  than  aligned,  then  the  maxi¬ 
mum  embedment  length  approaches  t/2  in.  Similarly,  the  n  m  refers  to 

max 

the  crack  surface  displacement  where  zero  aggregate  bridging  and  inter¬ 
locking  pressure  is  assumed  (Fig.  6). 

If  the  stresses  in  these  nonlinear  zones  (fiber  bridging  and  matrix 
process  zone)  are  assumed  to  be  purely  under  uniaxial  tensile  behavior, 
then  crack  length  "a"  can  be  replaced  by  an  effective  (elastic)  crack 

such  that  aefj  ■  a  +  *>p,  where  ^  is  the  idealized  length  of  the 
matrix  process  zone  (Fig.  5).  This  effective  crack  "ae£f"  sustains  two 
types  of  crack  closing  pressures:  one  due  to  the  fiber  bridging  and  another 
due  to  aggregate  bridging  which  was  neglected  as  explained  earlier.  The 
concept  of  this  model  is  somewhat  similar  to  that  originally  proposed  by 
Dugdale  [lO]  and  Barenblatt  [18] . 

For  a  given  applied  load  and  an  effective  crack  length,  the  crack  open¬ 
ing  process  is  primarily  resisted  by  the  bridging  of  fibers  across  the  crack. 
It  is  necessary  to  first  calculate  the  size  of  the  fiber  bridging  zone. 

Since  it  is  a  function  of  the  crack  profile  and  the  length  of  the  fibers 
used,  an  approximate  crack  profile  was  first  assumed  (a  linear  crack  pro¬ 
file  was  assumed  In  this  study).  Knowing  the  crack  profile,  crack  length 

and  n*  (half  the  fiber  length),  the  size  of  fiber  bridging  zone  was  then 
max 

determined.  From  the  uniaxial  tensile  stress-displacement  relationship  and 
the  assumed  crack  profile,  the  fiber  bridging  pressure  can  also  be  obtained. 
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The  stress  displacement  relationships  for  different  fiber  volume  fraction 
and  aspect  ratio  are  given  in  Table  1.  Details  of  how  to  obtain  such  re¬ 
lationships  will  be  discussed  later. 

In  this  analysis,  it  is  assumed  that  crack  will  Initiate  as  the  crack 

surface  displacement  at  the  tip  of  the  effective  crack  length  "a"  reaches 

a  value  (defined  in  Figures  5  and  6). 

max 

As  the  fiber  bridging  pressure  distribution  depends  on  the  crack  sur¬ 
face  displacement,  which  in  turn  is  a  function  of  the  applied  load,  speci¬ 
men  geometry,  the  size  of  fiber  bridging  zone,  the  length  of  the  matrix 
process  zone  A  and  the  closing  pressure  itself,  an  iterative  procedure 
was  then  needed  in  the  analysis  as  follows: 

Consider  a  given  crack  length  "a"  in  a  fiber  reinforced  concrete 
specimen  just  prior  to  its  further  extension. 

1.  Assume  a  crack  profile  and  a  matrix  process  zone  of  length  Ap. 

With  the  given  length  of  fibers  used,  the  size  of  fiber  bridging 

zone  A  can  be  calculated  using  1 1  ^  as  the  limit  of  the  fiber 
f  max 

bridging  zone. 

2.  Knowing  A^  and  the  assumed  crack  profile,  calculate  the  closing 
pressure  distribution,  using  the  stress-displacement  relation¬ 
ship  given  in  Table  1.  Note  that  for  simplicity  the  closing 
pressure  distribution  in  the  matrix  process  zone  is  neglected. 

3.  For  a  given  specimen  geometry  (double  cantilever  beam  and  notched 
beam  specimens  were  considered  here) ,  the  applied  load  P  and  the 
crack  closing  pressure,  calculate,  using  theory  of  elasticity,  the 
crack  opening  displacements  for  the  effective  crack  aeff. 


k 

.  a 

w 

J, 
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4.  If  the  crack  opening  displacement  at  the  end  of  matrix  process 

rone  (  £_)  Is  equal  to  n  m  ,  then  the  Initiation  criterion  is 
P  max 

satisfied  and  the  assumed  value  of  is  partially  a  correct  one; 

otherwise  a  new  value  of  £  is  assumed  and  the  above  steps  are 

P 

repeated  until  the  initiation  condition  is  satisfied. 

5.  To  further  ensure  that  the  iterated  £  and  the  assumed  crack 

P 

profile  are  correct,  the  load  line  deformation  (  P  )  is  cal¬ 
culated  based  on  the  same  elastic  principle  and  then  compared 
with  the  experimentally  observed  values.  If  these  values  do  not 
correspond  to  the  measured  data,  a  new  crack  profile  is  assumed 
and  the  above  procedure  is  repeated  until  this  condition  is 
satisfied. 

NORMALIZED  UNIAXIAL  TENSILE  STRESS-DISPLACEMENT  RELATIONSHIP 

Experiments  with  fiber  reinforced  concrete  specimens  subjected  to  uni¬ 
axial  tension  in  a  rigid  testing  system  have  shown  that  the  post-peak  dis¬ 
placements  of  the  specimens  are  essentially  a  result  of  the  opening  of  a 
single  crack  [22].  Thus  the  uniaxial  tensile  stress-displacement  relation¬ 
ship  in  the  post  cracking  region  can  be  taken  as  the  crack  closing  pressure 
(*)  versus  crack  surface  displacement  (n  )  relationship. 

Since  the  uniaxial  tensile  stress-displacement  is  important  to  the 
fracture  behavior  of  fiber  reinforced  concrete,  careful  study  of  this  rela¬ 
tion  is  essential,  particularly  in  the  post  cracking  region.  Many  inves¬ 
tigations  have  conducted  uniaxial  tensile  tests  of  fiber  reinforced  con¬ 
crete  with  different  fiber  volume  fractions  and  aspect  ratios  [l2,  22,  23J. 
These  results  were  used  for  the  analysis  here. 
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The  stress-displacement  relationship  of  fiber  reinforced  composites  in 
the  post-cracking  region  is  bounded  by  two  parameters  (Fig.  7):  the  maximum 
post-cracking  strength  (  and  the  maximum  pulled-out  displacement  of 

flber  <,iax)* 

The  maximum  post-cracking  strength  is  generally  known  to  be  a  function 
of  fiber  volume  fraction  (V^),  aspect  ratio  (———),  fiber  distribution  and 
embedment  length  coefficient  (  0  or  4  in  Eq.  4-a)  and  the  bond 

strength  of  fiber  matrix  interface  ( r  )  £l5]  where: 

-  Q  r  V.  (15) 

"max  M  f 

The  maximum  fiber  pulled-out  displacement  normally  depends  on  the  em¬ 
bedment  length  and  the  fibers  orientation.  It  can  be  seen  (Fig.  6)  that 
fibers  may  bridge  over  a  crack  in  different  ways;  some  fibers  may  have  most 
of  the  length  embedded  on  one  side,  and  have  only  a  small  portion  left  on 
the  other  side  of  the  crack  (labelled  1);  some  fibers  may  be  totally  em¬ 
bedded  on  one  side  of  the  crack  (labelled  3) ;  and  hardly  but  possibly  that 
a  fiber  may  be  equally  embedded  on  both  sides  of  the  crack  (labelled  2  in 
Fig.  6).  Since  the  smaller  side  of  the  embedment  length  is  the  one  which 
will  be  pulled  out  of  the  matrix  (this  is  due  to  a  smaller  total  resisting 
force  or  frictional  force),  it  can  then  be  concluded  that  the  maximum 
pulled-out  displacement  cannot  exceed  half  the  fiber  leneth  which  is  used 

here  as  the  value  of  n  . 

max 

By  normalizing  the  post-cracking  stress  with  the  maximum  post-cracking 

strength  (  */*aax)  and  the  post-cracking  displacement  with  (~~ — ■)  , 

’’max 

it  was  found  that  there  exists  a  unique  relationship  between  the 
normalized  post-cracking  stress  and  displacement  for  a  given  type  of  fiber. 
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This  relationship  was  independent  of  the  fiber  volume  fraction  and  aspect 
ratio  as  can  be  seen  in  Fig.  8.  To  generalize  a  unique  normalized  stress- 
displacement  relationship  for  steel  fiber  reinforced  concrete,  it  was 
assumed  that  the  maximum  post-cracking  strength  occurs  at  the  end  of  the 
matrix  process  zone  or  the  beginning  of  the  fiber  bridging  zone,  and 
smoothly  decreases  to  zero  as  the  crack-surfaces  displacement  reaches  half 
of  the  fiber  length.  Therefore,  the  boundaries  for  this  relationship  can 


be  summarized  as: 

at  x  *  o  or  q  ■  o  :  — — —  *  1.0 

*max 

at  x  -  *■,  or  n  ■  1  ^  :  — 2 —  «  0.0 

£  max 


(16) 


where  x  is  the  distance  from  the  crack  tip,  n  is  the  crack  surface  dis¬ 
placement  at  x.  For  the  condition  of  monotorically  decreasing  stress,  we 
can  write: 

-4~—  <  o  for  o  <  i|  <  i|^  (17) 

u  n  max 

or  o  <  x  < 

At  attempt  to  predict  the  normalized  stress-displacement  relationships 
for  steel  fiber  reinforced  concrete  based  on  the  experimental  results  re¬ 
ported  in  [ll,  22,  23]  is  made  here.  To  satisfy  all  the  required  boundary 
conditions  given  in  Eqs.  (16,  17),  one  of  the  possible  solutions  is: 


where  9  is  the  maximum  post-cracking  strength  and  q  is  the  maximum 
max  max 

pulied-out  displacement  of  fibers. 


For  ocher  Cypes  of  fiber  reinforced  composites,  Che  normalized  stress- 
displacement  relationship  given  in  Eq .  18  may  be  modified  with  an  exponen¬ 
tial  factor  of  — l]B  where  Eq.  18  can  be  rewritten  as: 


where  m  and  n  are  constants  which  will  depend  on  the  type  and  the  pull-out 
behavior  of  fibers. 

It  can  be  seen  that  the  post-cracking  strength  at  any  pulled-out  dis¬ 
placement  ( V  )  can  be  predicted  provided  that  p max  and  p  are  known 
(  is  a  function  of  P^  ,  Pg  ,  T  ,  Vf ,  and  i/d  while  p  equals 

half  the  fiber  length).  Thus  the  normalized  relationship  shown  in  Eqs.  8, 
18,  19  could  be  very  useful  in  design  of  fiber  reinforced  concrete  com¬ 
posites  . 

To  predict  the  tensile  stress-displacement  relationships  in  this 
study,  the  value  of  A^  in  Eq.  4-b  was  neglected  and  the  factor  A2  which  is 
a  function  of  ^  ,  ^g  ,  and  r  has  been  selected  from  the  literatures 
[15,  24  -  26]  as  follows: 

For  the  double  cantilever  beam  specimen,  with  the  method  of  casting, 
the  type  of  fibers  used  and  the  critical  specimen  thickness  of  0.5  in.  M. 
the  value  of  A2  was  selected  as  600  psl.  In  the  case  of  the  notched  beam 
specimens,  the  different  fibers  size,  the  method  of  casting,  and  specimen 
thickness  of  1.5  in.  [6],  lead  to  the  value  of  A2  as  300  psl.  With  these 
selected  values  of  A2,  the  stress-displacement  relationship  can  be  predic¬ 
ted,  using  Eqs.  15  and  18. 
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Different  a~n  relationships  for  different  fiber  volume  fraction 
and  aspect  ratio  used  in  this  study  are  given  in  Table  1. 

CALCULATION  OF  CRACK  SURFACE  DISPLACEMENT 
Most  fracture  studies  for  cementitious  composites  in  the  past  were 
based  on  the  calculation  of  the  stress  intensity  factor  K  using  the  linear 
elastic  fracture  mechanics  concept.  In  this  study,  a  method  is  proposed 
to  calculate  the  strain  energy  release  rate  G  based  on  the  change  of  global 
compliance  — — —  approach. 

Approximate  and  simple  methods  of  calculating  compliances  of  a 
cracked  double  cantilever  beam  and  centrally  notched  beam  are  described  in 
the  following  section. 


DOUBLE  CANTILEVER  BEAMS 


The  compliance  of  the  double  cantilever  beam  specimens  was  calculated 
using  the  so-called  DCB  approach  [19,  20j. 

A  single  cantilever  beam  model  under  the  applied  load  P  and  the  crack 
closing  pressure  is  shown  in  Fig.  9.  The  deformation  for  the  double  can¬ 
tilever  beam  in  this  analysis  was  calculated  from  the  beam  bending  theory 
where  deformation  due  to  the  applied  point-load  P  for  a  single  cantilever 
beam  (Fig.  9)  can  be  determined  from  [l9]: 


ypM 


12P 

BE 


aeffx 


1  .  H(x) 

-  -j  m - <- 


2H  (x)H2  (o)  C2  H<°> 


1  x 

2  2 


+  1 


l  i 

+  4  7T 


in  MiS*. 
H(x) 


H(x)  H(o) 


Co  ■  0 


(20) 


where  yp(x)  is  the  deformation  at  any  point  x  from  the  effective  crack  tip, 


aeff*  c’  and  C2  are  Che  beam  depths  at  the  free  end  and  the  slope  of  the 
tapered  beam.  E  and  u  are  material  constants,  and  B,  H(x)  and  H(0)  are 
beam  width  and  thicknesses  at  distance  x  from  crack  tip  and  at  the  crack 
tip.  These  variables  are  shown  in  Fig.  9. 

For  the  deformation  due  to  fiber  bridging  pressure,  with  the  nonlinear 
crack  closing  pressure,  a  numerical  integration  was  employed  to  solve  the 
beam  bending  governing  equation: 


d2y 

~dxr 


M 

El 


AG 


dV 

dx 


(21) 


where  A  and  I  are  cross  sectional  area  and  moment  of  inertia,  E  and  G  are 
elastic  and  shear  modulus  while  K  is  the  average  shear  constant.  M.ahd  V 
are  moment  and  shear  which  are  given  as: 

for  ^  <  x  <  ^  ; 

P  f  P  t  l 

f  +  p 

M(x)  _  i 

tn.  <r( t)  (t-x)  dt 


V(x) 


l 

l 

r 


(22) 


*(t)  dt 


for  o  <  x  <4  . 

P  * 


l  l 


M(x) 


-n 


/^f  +  ^p 
c 


V(x) 


for  x  >£f  +  4p  . 


»(t)  (t-x)  dt 


»(t)  dt 


(23) 


M(x)  -  0  ;  V(x)  -  0 


(24) 
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The  total  deformation  is  equivalent  to  the  deformation  due  to  annlied 
load  P  subtracted  by  those  due  to  fiber  bridging. 

(np^ total  "  (’Mapolied  load  ~  < n P> fiber  bridging  (25) 

Details  of  the  test  setun  and  snecimen  dimensions  are  given  in  Fig.  10. 


CALCULATION  OF  CRACK  MOUTH  DISPLACEMENT  FOR  XB  SPECIMENS 


A  fiber  reinforced  concrete  beam  with  an  initial  cast  notch  of  length 

•  ’ 

a0,  subjected  to  pure  bending,  is  shown  in  Fig.  11.  As  similar  to  the  case 
of  double  cantilever  beam  analysis  for  a  given  load  P  and  crack  length  *a", 
all  nonlinear  zones  (matrix  process  zone  and  fiber  bridging  zone)  were  re¬ 
placed  by  a  region  of  crack  closing  pressure.  The  crack  closing  pressure 
within  the  matrix  process  zone  was  assumed  to  be  small  compared  to  those 
due  to  fiber  bridging  and  therefore  was  neglected. 

The  fiber  bridging  stress  again  defends  on  the  assumed  crack  profile 
which  is  a  function  of  crack  growth  "a"  and  crack  mouth  displacement  "CMD". 
Using  the  normalized  stress-displacement  relationship  given  in  Eq.  18,  the 
calculated  maximum  post-cracking  strength  (  *max)  anc*  the  assumed  crack 
profile,  the  fiber  bridging  stress  distribution  can  then  be  determined  as 
follows : 

For  the  assumed  linear  crack  profile: 


CMD 


and 


_n_ 

y 

CMD 


Note  that  Eq.  26  gives  n  “  0  when  v  ■  0  which  is  in  fact  not  true, 
for  the  simplicity  in  calculating  the  fiber  bridging  pressure,  this 


(26) 

Only 
as sump - 


•  .*  -• V'  >  - 


tion  is  made. 
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Substitute  H  from  Eq.  26  into  the  normalized  stress-displacement  re¬ 
lationship  (Eq.  19).  Eq .  18  can  be  rewritten  as: 


a  (y)  = 


’  1  -  — i- 

jx  f 

n  n< 


y 

a 


where 


is  the  maximum  oulled-out  disolacement  (equal  to  half  the  fiber 


length) . 


By  introducing  »  ( - j - )  ,  Eq.  27  can  be  rewritten  as: 

a  •  H  max 


I1 ' v] 


Knowing  the  fiber  bridging  stress  acting  along  the  length  of  on 
the  effective  crack  of  length  a  +  ,  the  change  in  rotation  of  the  beam 

resulting  from  the  effective  crack  -  aeff  can  be  calculated  from  the  rela¬ 
tions  derived  by  Okamura  et  al  [27,  28]  : 

9  -  X^  .  F  +  XM  .  M  (29) 


where  9  is  the  relative  rotation  of  the  crack  front,  F  is  the  axial  force 

resulted  from  the  fiber  bridging  pressure,  M  is  the  resultant  moment  due  to 

the  external  load  and  those  resisting  moment  due  to  F  (Fig.  11),  and  X^ 

and  X^  are  the  increments  of  the  compliances  caused  by  the  presence  of  the 

effective  crack  -  a  and  are  given  as: 
err 


X -  2(1-  y2)  _6_  f 
™  BE  W  J 


t  Yy  (  S  )  Ym  (  i  )  d  £ 


(30-a) 


2(1-  \ 
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*  *■/■ 


{  ym  (  O  d  { 
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where  5  ■  -  and  ^f($)»  V«(  £  )  are  given  as: 

W 

Y_( £  )  -  1.99  -  0.41(  £  )  +  18.70( ?  )2  -  38.48( £  )3  +  53 . 85 (  £  )4 

P 


(30 -c) 


Y„(  £  )  -  1.99  -  2.47(  £  )  +  12.97(  £  )2  -  23.17( £  )3  +  24.80(  £  )4 

M 


(30-d) 


Note  that  Eqs.  29  and  30  are  applicable  only  for  £  S  0.6. 


For  the  case  of  non-uniformly  distributed  fiber  bridging  stress,  the 
total  F  can  be  determined  from: 


F  -  B 


a  (y)  dy 


substitute  o(y)  Eroin  E9*  28 


F  ■  B  a 


F  "  B  °max  • 


/  '  t1  ‘  VI 

*  l{  \  li 


2  t  1 

\  lf +  \  — — 

3 


Note  that- the  size  of  fiber  bridging  zone  was  again  controlled  by  half  the 
fiber  length;  however,  in  the  case  of  notched  beam  specimen  the  value  of 
COD  was  always  smaller  than  half  the  fiber  length.  As  a  result,  the  size 
of  &£  was  equal  to: 

for  COD  ^  L/2  I 

tf  -  a  -  a0  (34) 

The  resisting  moment  due  to  force  F  depends  on  the  distance  between 
the  neutral  axis  and  the  crack  tip.  The  problem  then  depends  on  where  the 
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location  of  the  neutral  axis  is.  According  to  Hannant  [l5]  ,  the  neutral 

axis  of  an  unnotched  fiber  reinforced  concrete  beam  is  approximately  one- 

ouarter  of  the  beam  depth  (or  the  unnotched  distance)  from  the  compressive 

surface.  Therefore,  the  distance  between  the  crack  tip  and  the  neutral  axis 

was  assumed  to  be  equal  to  (_ W  +  — - a) 

4  4  ' 


The  resisting  moment  can  then  be  determined  from: 


Mr  "  Bfl  a(y)  |^(-|-  W  +  ~~  *  «)  +  y  J 


3  a 

Introducing  -  ~£~  W  +  °  -  a  ,  and  substituting  a(y)  from  Eq.28: 


r 

”t  "  BJ  0 »  U  -  V>  <dm  +  y)  dy 
0 

4-  if2  (  A  2  j  o  .  \  2.3 


+  <  *k  «m  -  2  \  )  + 


Ak 


For  a  given  value  of  the  applied  bending  moment,  the  dimensions  of  the 

beam  and  the  length  of  the  notch,  a  value  of  can  be  iterated  such  that 

the  calculated  value  of  0  from  Eq .  29  is  equal  to  H  ”  /  *■_  (the  value  of 

max  h 

_2 

ti  m  was  considered  from  the  matrix  and  was  found  to  be  0.08  x  10  in. 

max 

according  to  [l9,  20,  2l]  ). 

Once  the  initiation  condition  is  satisfied,  the  crack  mouth  displace¬ 
ment  is  calculated  from: 

CMD  -  0  (a  +  *■  )  (38) 

P 
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the  calculated  value  of  CMD  is  then  checked  with  the  experimentally  ob¬ 
served  data.  If  this  condition  is  satisfied,  the  assumed  crack  profile 
and  the  iterated  are  correct,  otherwise  a  new  crack  profile  is  assumed 
and  the  above  procedure  is  repeated  until  such  displacement  condition  is 
satisfied. 

Details  of  the  test  setup  and  the  specimen  dimensions  for  the  notched 
beam  specimen  are  given  in  Fig.  12. 

CALCULATION  OF  FRACTURE  ENERGY 

The  resistance  to  crack  extension  can  be  expressed  in  terms  of  the 
strain  energy  release  rate.  This  term  should  include,  in  addition  to  the 
elastic  energy,  the  energv  absorbed  in  the  fiber  bridging  zone  (as  well  as 
in  the  matrix  process  zone,  which  was  ignored  here).  In  terms  of  the 
global  compliance  concept,  the  terms  for  strain  energy  release  rate  should 
include  both  the  unloading  compliance  (C^  in  Fig.  13)  and  the  permanent 
deformation  (  5p) . 

To  determine  the  actual  amount  of  energy  required  for  crack  propaga¬ 
tion,  it  has  been  shown  that  the  classical  concept  of  strain  energy  release 
rate  must  be  modified  to  include  the  effect  of  nonlinear  zone  (in  this  case, 
the  fiber  bridging  zone)  [7,  9,  19].  It  was  also  shown  that  the  energy 
required  for  incremental  crack  growth  should  be  the  area  under  the  load 
deformation  curve  between  two  unloading  lines  at  those  corresponding  crack 
growths  (shaded  area  in  Fig.  13),  i.e.. 


(39) 
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where  Pj  and  ?2  are  two  consecutive  neighboring  loads,  tn  is  the  critical 
thickness  in  the  crack  plane,  CR  and  ip  are  the  unloading  (elastic)  com¬ 
pliance  and  the  oermanent  deformation,  respectively.  To  calculate  Gg,  from 
Eq.  39,  it  is  clearly  indicated  that  the  change  of  unloading  compliance  and 
permanent  deformation  with  crack  growth  must  be  determined. 

During  the  unloading  process,  the  fiber  slipping  distance  or  pulled- 
out  displacement  was  assumed  to  remain  constant.  This  behavior  was  exper¬ 
imentally  observed  by  Hawkins,  Lin  and  Jeang  [29].  They  reported  that  the 
slip  in  reinforced  concrete  during  unloading  remains  constant  (Fig.  14). 

It  Implies  that  the  same  fiber  bridging  pressure  which  resists  the  crack 
from  opening  during  the  loading  process  also  acts  as  frictional  force  and 
resists  the  crack  from  closing  during  the  unloading  process.  As  a  result, 
we  can  conclude  that  the  existing  permanent  deformation  at  the  end  of  the 
unloading  period  (P  ■  0)  is  equal  to  the  load-line  deformation  due  to  fiber 
bridging  stress  which  prevents  the  crack  from  opening.  If  the  fiber  bridg¬ 
ing  stress  was  assumed  to  remain  constant  throughout  the  unloading,  the 
deformation  2  (due  to  fiber  bridging  during  loading)  must  equal  the  defor¬ 
mation  3  (due  to  fiber  preventing  crack  from  closing  during  unloading) 

(Fig.  15). 

To  determine  R-curve  (fracture  energy  vs.  crack  extension)  from  a  given 
set  of  load  P  and  crack  a,  the  procedure  mentioned  earlier  for  DCB  and  NB 
was  used  to  Iterate  for  the  size  of  the  matrix  process  zone  and  calculate 
the  total  load-line  deformation  (^p  )  total.  From  the  same  crack  closing 

pressure,  the  permanent  deformation  (Hp  )  perm  can  also  be  calculated  from 
the  elastic  beam  bending  theory.  Knowing  both  the  oermanent  deformation 
and  the  total  load-line  displacement,  the  elastic  deformation  can  be  deter¬ 
mined  from  Fig.  15: 


-  --- 
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C  Tip)  total  -  (  Up)  perm.  +  (  Tip)  elastic  (40) 

The  unloading  compliance  (CR)  can  then  be  calculated  from: 

CR  -  <  Viatic  (41) 

By  clotting  C  and  ( h p)  perm,  with  crack  length,  the  change  of  un¬ 
loading  compliance  and  permanent  deformation  with  respect  to  crack  growth 
can  be  determined.  Consequently,  Gg  can  be  calculated  from  Eq.  39,  and 
eventually  R-curve  can  be  obtained  by  plotting  the  obtained  Gp  with  the 
actual  crack  growth, 


EXPERIMENTAL  INVESTIGATION 


SPECIMEN  CONFIGURATION 

The  double  cantilever  beam  specimens  (Fig.  10)  used  in  this  study 
were  24  in.  long  tapered  beams  developed  by  Visalvanich  and  Naaman  [ll] . 
The  specimens  were  2  in.  thick  with  double  grooves  along  the  center  line 
of  the  specimen  which  provided  a  0.5  in.  critical  thickness.  Details  of 
the  specimen  setup  were  given  in  [11]. 

The  notched  beam  specimens  were  tested  by  Velazco  et  al  (30) .  The 
specimens  were  1.5  in.  thick  with  3  in.  width  and  a  span  length  of  15  in. 
Details  of  the  specimen  and  test  setup  are  shown  in  Fig.  12. 

MATERIAL  PROPERTIES  AND  MIX  PROPORTION 

The  matrix  mix-proportion  used  in  this  study  was  (1:2:0. 5)  (Cement  : 
Sand  :  Water) .  The  cement  used  was  high  early  strength  ASTM  Type  III  ce¬ 
ment;  the  sand  was  a  siliceous  sand  with  a  maximum  particle  size  passing 
sieve  No.  8  (0.0937  in.  opening). 

The  fibers  used  were  straight  cut,  brass  coated  steel  fibers  with  the 

3 

specific  gravity  of  490  lb/ft  .  Three  different  volume  fraction  of  fibers, 
0.5,  1  and  2  percent  were  selected  to  study  the  effect  of  fiber  content. 

Two  different  sizes  of  fiber:  0.25  in.  long  with  0.006  in.  diameter,  and 
0.75  in.  long  with  0.016  in.  diameter,  were  used.  In  this  report,  four 
series  of  fiber  volume  fraction  and  'aspect  ratio  (as  listed  in  Table  1) , 
with  two  specimens  each,  were  theoretically  studied.  Only  the  series  of 
Vf  »1Z  with  1  ■  0.25",  and  d  ■  0.006"  which  were  tested  by  the  authors 
with  the  loading-unloading  technique  [7,  19];  all  other  series  of  DCB 
specimens  were  tested  by  Visalvanich  [ll]  without  the  unloading  process 
and  therefore  could  not  be  used  to  obtain  the  experimental  R-curves. 
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The  notched  beam  tests  were  conducted  by  Velazco  [30] .  The  data  used 

in  this  study  have  a  constant  fiber  volume  fraction  of  1  percent  with  the 

Initial  cast  notch  varied  from  3/8  in.  to  1.5  in.  The  steel  fibers  used 

were  brass-coated,  1  in.  long,  and  0.01  in.  in  diameter  (  m  100).  The 

■» 

fibers  specific  gravity  was  490  lb/ft  ,  and  the  matrix  mix-proportion  was 
1:2:0. 5  (Cement  :  Sand  :  Water).  The  experimental  data  provided  by  Velazco 
[30]  were  the  applied  load,  crack  growth  and  the  crack  mouth  displacement; 
no  unloading  was  reported  in  the  testing  procedure. 
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COMPARISON  WITH  EXPERIMENTAL  DCS  DATA 

Since  only  the  series  of  V  ■  1  %,  i  *  0.25",  and  d  *  0.006"  for  the 
double  cantilever  beam  specimens  were  tested  with  the  loading-unloading 
technique,  only  the  experimental  R-curves  of  these  series  were  then  cal¬ 
culated  and  comnared  with  the  theoretical  values  in  this  section,  while 
the  results  of  all  other  series  will  be  reported  later. 

With  the  assumed  crack  profile,  the  iterated  Ep  was  calculated  by 

satisfying  the  initiation  criterion  (i.e.,  crack  displacement  at  the  end  of 

the  matrix  process  zone  enuals  0  m  of  the  matrix,  Figs.  5  and  6).  To 

max 

ensure  that  the  obtained  and  the  assumed  crack  profile  were  correct, 
the  predicted  total  load-line  deformation  must  be  in  good  agreement  with 
the  observed  values. 

Fig.  16  shows  the  plot  of  total  load-line  deformation  (np)  with 

crack  growth  while  Fig.  17  reports  the  load-deformation  curves.  It  can  be 
seen  that  the  predicted  values  are  in  good  agreement  with  the  experimentally 
observed  values.  (Note  that  the  solid  points  refer  to  the  experimental  data 
while  the  blank  point  represents  the  predicted  values;  the  same  notation 
will  be  used  throughout  this  study.) 

Since  the  fiber  bridging  pressure  primarily  resulted  from  the  fric¬ 
tional  bond  strength  and  the  slipping  distance  was  assumed  to  remain  con¬ 
stant  during  unloading  [29] ,  the  permanent  deformations  can  easily  be  deter¬ 
mined  from  the  obtained  fiber  bridging  pressure,  using  the  elastic  beam 
bending  concept.  Fig.  18  shows  a  clot  of  permanent  deformation  versus 
crack  growth  and  it  can  be  seen  that  the  predicted  results  are  in  good 
correlation  with  the  experimental  data. 
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In  Fig.  19,  a  good  correlation  is  also  observed  for  the  predicted  and 
measured  unloading  compliances  at  small  crack  growth.  The  correlation  is 
poorer  at  the  larger  crack  growth.  This  may  be  because  of  the  so-called 
end-effect.  The  end-effect  in  the  double  cantilever  beam  specimen  is  re¬ 
lated  to  the  size  of  the  uncracked  portion  of  the  specimen  since  the  DCB 
analyses  assume  the  beam  to  be  fixed  at  the  crack  tip.  As  a  result,  the 
smaller  the  uncracked  ligament,  the  more  inaccurate  the  fixed-end  concept. 
For  the  double  cantilever  beam  used  in  this  study,  the  valid  region  without 
the  end-effect  is  for  crack  growth  of  less  than  16  in. 

Knowing  the  change  of  permanent  deformation  and  unloading  compliance 
with  crack  growth  (Figs.  18  and  19),  R-curves  can  be  determined  from  Eq. 

1 391  The  predicted  R-curve  was  compared  with  the  exnerlmentally  observed 
R-curve  in  Fig.  20  where  good  agreement  was  observed  in  general. 

COMPARISON  WITH  OTHER  DCB  DATA 

Other  DCB  experimental  data  which  will  be  used  for  comparison  with  the 
predicted  results  are  the  series  of  Vf  -  0.5  X,  1  -  0.25",  d  -  0.006"; 

Vf  -  2  X,  i  -  0.25",  d  -  0.006";  and  Vf  -  1  X,  *  -  0.75",  d  -  0.016". 

The  comparisons  of  the  theoretically  predicted  load-line  deformation 
with  the  experimental  results  reported  in  [11]  are  presented  in  Figs. (21a  - 
21c)  which  are  the  plot  of  load-line  deformation  with  crack  growth.  It  can 
be  seen  that  all  predicted  results  were  in  good  correlation  with  the  exper¬ 
imental  data.  This  implies  that  the  assumotlons  of  linear  crack  profile  and 
the  negligibility  of  crack  closing  pressure  within  the  matrix  process  zone 
are  satisfactory.  It  was  also  found  that  the  higher  the  fiber  volume  frac¬ 
tion  (Vf ) ,  the  larger  the  total  load-line  deformation. 
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The  load-deformation  carves  predicted,  using  the  proposed  model,  were 
found  to  be  in  good  agreement  with  the  experimentally  observed  values  re¬ 
ported  in  [ll] ,  as  can  be  seen  in  Fig*.  22a  -  22c. 

Fig.  23  shows  a  plot  of  predicted  permanent  deformation  for  all  four 
sets  of  fiber  reinforced  concrete.  It  was  observed  that  the  higher  the 
volume  fraction,  the  longer  the  predicted  permanent  deformation.  It  can 
also  be  seen  that  the  larger  the  value  of  fiber  aspect  ratio,  the  larger 
the  permanent  deformation.  The  rate  of  change  in  permanent  deformation 
(with  respect  to  crack  growth  also  increased  with  fiber  volume  fraction 
and  aspect  ratio. 

A  plot  of  unloading  compliance  with  crack  growth  for  all  four  sets  of 
fiber  reinforced  concrete  was  shown  in  Fig.  24.  For  Vj  *  0.5  %,  the  pre¬ 
dicted  C^-a  relationship  was  close  to  linear.  It  was  also  observed  that 
the  larger  the  amount  of  fibers  (higher  volume  fraction) ,  the  more  non¬ 
linear  behavior  exists  for  the  unloading  compliance  and  crack  growth  rela¬ 
tionship  . 

Fig.  25  shows  the  comparison  of  the  predicted  R-curves  for  all  four 
series  of  fiber  reinforced  concrete.  It  can  be  seen  that  increasing  the 
amount  of  fibers  (volume  fraction)  resulted  in  an  increase  in  the  fracture 
resistance.  For  the  same  volume  fraction,  increasing  the  fiber  aspect 
ratio  (length)  also  increased  the  fracture  toughness  of  the  composites. 

Note  that  if  the  asymptotic  value  of  the  R-curve  is  considered  as  a 
material  parameter,  then  that  value  can  be  a  useful  quantity  in  identifying 
the  benefit  of  fiber  addition.  For  example,  by  comparing  the  R-curves  of 
specimens  reinforced  with  0.75  in.  long  fibers  and  plain  mortar,  it  can  be 
observed  that  the  asymptotic  steady-state  value  of  fracture  energy  is 
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Increased  more  than  40  times  due  to  fiber  addition.  This  relative  value  is 
somewhat  comparable  with  the  reported  value  of  "toughness  index"  (relative 
values  of  areas  under  the  load-deflection  curve  in  flexure)  by  other  inves¬ 
tigators  £33,  34,  353.  Thus,  the  model  pronosed  here  seems  to  provide  a 
tool  for  predicting  the  effectiveness  of  fibers  in  imoroving  the  fracture 
resistance  of  concrete. 

COMPARISON  WITH  EXPERIMENTAL  NOTCHED  BEAM  DATA 

• 

Using  the  theoretical  model  proposed  for  the  notched  beam  specimens, 
the  theoretical  load-crack  mouth  displacement  relationships  were  predicted 
for  nlain  mortar  and  fiber  reinforced  concrete  (Fig.  26).  The  mix-propor¬ 
tion  for  plain  mortar  was  1:2:0. 5  (Cement  :  Sand  :  Water).  Three  different 
sets  of  initial  notch  aQ  with  Vf  equalled  to  1  percent  and  "amax  -  300  psi" 
were  studied. 

A  good  correlation  is  observed  for  the  predicted  and  the  measured 
crack  mouth  displacement  for  plain  mortar  (Fig.  26).  These  results  were 
calculated  using  the  observed  peak  loads  and  the  initial  crack  length  from 
[30]. 

For  fiber  reinforced  concrete,  the  comparisons  in  Fig.  26  show  that 
the  analysis  does  not  correlate  exactly  with  the  experimental  results. 
However,  the  predicted  values  were  still  in  good  agreement  with  the  ob¬ 
served  data.  These  variations  may  be  attributed  to  the  fact  that:  1)  the 
actual  neutral  axis  for  the  notched-fiber  reinforced  beam  is  not  known,  the 
neutral  axis  used  in  this  analysis  was  based  on  what  has  been  proposed  by 
Hannant  [15]--  (d/^  from  the  compression  face  where  d  is  the  depth  of  the 
unnotched  beam);  2)  the  theoretical  concept  proposed  by  Okamura  et  al 
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£27,  28]  is  valid  only  when  the  a/w  ratio  is  less  than  0.6,  with  fiber  re¬ 
inforced  concrete,  most  observed  crack  growths  reoorted  in  Ref.  [30]  had 


a/w  much  larger  than  0.6. 


The  assumed  linear  crack  profile  might  not 


be  accurate  enough  for  the  case  of  notched  beam  specimens  which  leads  to  a 
wrong  estimation  of  the  fiber  bridging  pressure. 

Additional  theoretical  analysis  for  predicting  crack  mouth  displace¬ 
ment  more  accurately  for  notched  beam  specimens  is  currently  under  way. 


MATRIX  PROCESS  ZONE.  FIBER  BRIDGING  ZONE  AND  CRITICAL  COD 


The  value  of  the  matrix  process  zone  4^  can  only  be  obtained  through 
the  iteration  procedure  which  must  satisfy  the  initiation  criterion  that 
the  crack  tip  opening  displacement  equals  H  ™  used  here  for  the  matrix 
was  0.08  x  lO-^  in.,  as  reoorted  in  Ref erences [l9 ,  20,  and  2l],  The  fiber 
bridging  zone  was  calculated  from  the  assumed  crack  profile  with  the  condi¬ 
tion  that  fibers  will  be  completely  pulled  out  when  the  crack  surface  dis¬ 
placement  equals  r|^max  (half  the  fiber  length,  Fig.  6).  In  this  study, 
since  the  value  of  crack  opening  displacement  was  always  less  than  r|fmax 
(half  the  fiber  length) ,  the  size  of  the  fiber  bridging  zone  was  then  equal 
to  the  crack  length  subtracted  by  the  length  of  the  initial  cast  notch. 

Fig.  27  shows  the  plot  of  both  the  calculated  matrix  process  zone  and 
the  fiber  bridging  zone  (  4p  and  4^,  respectively)  for  the  DCB  specimens. 

It  can  be  observed  that  the  matrix  process  zone  is  essentially  constant  with 
respect  to  crack  growth.  The  value  of  4^  calculated  from  fiber  reinforced 
mortar  specimen  was  almost  identical  to  those  obtained  from  mortar  specimens 
[19,  20,  2l].  The  conclusion  of  constant  process  zone  was  also  confirmed  in 
other  materials  like  DMMA  [31,  32].  The  obtained  value  of  4  was  found  to 


be  approximately  about  3  in.  and  seemed  to  be  independent  of  different  fiber 
volume  fraction  and  aspect  ratio  (Fig.  27). 

A  plot  of  the  matrix  process  zone  and  the  fiber  bridging  zone  with  the 
uncracked  ligament  in  the  notched  beam  snecimens  is  shown  in  Fig.  28.  The 
predicted  matrix  process  zones,  using  the  proposed  model,  were  found  close 
to  what  has  been  reported  for  the  unreinforced  mortar  [19  ,  20].  Fig.  18 
also  shows  that  the  fiber  bridging  zone  increases  with  crack  growth,  or  as 
the  uncracked  ligament  decreases. 


Table  1. 


PREDICTED  STRESS-DISPLACEMENT  RELATIONSHIPS 


FOR  DIFFERENT  Vf  AND  *7d  FOR  DCB  SPECIMENS 


vf 

* 

l 

(in.) 

d 

(in. ) 

a  _  p  relationship  ** 

(nsi) 

1 

0.5 

41.67 

0.25 

0.006 

<7  -  125  (1  -  80  )2 

1.0 

41.67 

0.25 

0.006 

o  -  250  (1  -  80  )2 

1.0 

46.88 

0.75 

0.016 

O  -  281.28  (1  -  2 .67  p  )2 

2.0 

41.67 

0.25 

0.006 

a  •  500  (1  -  8 n  ) 2 

These  stress-displacement  relationships  were  predicted,  using 
j}  r  *  600  psi,  and  the  proposed  normalized  relationship  (Eq.  15) 
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APPENDIX  I  -  DOTATIONS 

cross  sectional  area 

a  factor  for  composite  strength  that  related  to  the  extent 
of  matrix  cracking 

a  factor  for  composite  strength  that  related  to  fiber 
bridging  (  •  4  Hg  T  ) 

f 

a  constant  obtained  from  the  ratio  of  CMD/ (a  -  n  ) 

max 

crack  length 

traction-free  crack  length 

effective  crack  length  (  *  a  +  l  ) 

P 

beam  thickness 

boundaries  of  the  zone  of  uniform  crack  closing  pressure 

compliance  (  *  ^ /  ) 

P 

DCS  -  beam  depth  at  free  end 
slope  of  the  tapered  DCB  specimen 
unloading  -  reloading  compliance 
crack  mouth  displacement 
crack  opening  displacement 
diameter  of  fiber 

distance  from  crack  tip  to  the  assumed  neutral  axis 
modulus  of  elasticity 

total  axial  force  in  notched  beam  specimen 
shear  modulus  of  elasticity 
critical  strain  energy  release  rate 
modified  strain  energy  release  rate 
strain  energy  release  rate  -  resistance  curve 
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DCB  -  beam  depth  at  crack  tip 
DCB  -  beam  depth  at  x 
moment  of  inertia 
average  shear  constant 

stress  intensity  factor  due  to  applied  force  P 
critical  stress  intensity  factor 

stress  intensity  factor  at  crack  tip  calculated  without 
considering  the  presence  of  fibers 
stress  intensity  factor  for  crack  norpagation  in  fiber 
reinforced  matrix 

stress  intensity  factor  due  to  fiber  bridging 
stress  Intensity  factor  due  to  wires  bridging 
critical  stress  intensity  factor  under  mode  I  opening 
fiber  length 
zone  of  fiber  bridging 

mlcrocracked  zone  or  matrix  process  zone 
fiber  aspect  ratio 
applied  moment 

exponential  constants  for  variation  of  fiber  pulled-out 
behavior 
applied  load 

two  consecutive  neighboring  loads  in  the  loading  sequence 
critical  beam  thickness  in  the  crack  plane 
crack-surface  displacement  of  a  crack  without  fiber 


bridging 


ur  crack-surface  displacement  due  to  fiber  bridging 

V  shear  force 

Vj  volume  fraction  of  fibers 

Vm  percent  volume  of  matrix 

W  depth  of  notched  beam  specimen 

x  distance  from  crack  tip 

y  distance  from  crack  tip  in  fiber  bridging  zone  of  notched 

beam  specimen 

yp(x)  crack-surface  displacement  at  distance  x  from  crack  tip 
due  to  applied  force 

Yp( £ ) ,  Y^(  £  )  polynomial  functions  due  to  force  and  moment  in 
terms  of  £  for  notched  beam  specimen 
a  a  factor  (less  than  unity)  for  the  effective  process  zone 

8  efficiency  factor  due  to  fiber  distribution  and  embedment 

length 

o  uniform  crack  closing  pressure 

<Jm  tensile  strength  of  matrix 

°max  maximum  post-cracking  strength  of  fiber  reinforced  composites 
a(t)  non-uniform  crack  closing  pressure 

a(y)  crack  closing  pressure  at  distance  y  in  notched  beam  specimen 

(  a  )  maximum  strength  of  fiber  reinforced  composites 

c  max 

n  crack-surface  displacement 

efficiency  factor  for  the  embedment  length  of  fibers 

nQ  coefficient  for  fiber  distribution  and  orientation 

n  load-line  deformation 

P 

(  n  )  ,  .  elastic  load-line  deformation  (reversible) 

p  elastic 
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(  n  ) 
p  pern 


permanent  load-line  deformation  (irreversible) 


(  np)totai  total  load-line  deformation 

maximum  pulled-out  displacement  of  fiber  reinforced 

max 

matrix  (  ■  half  the  fiber  length) 


nm 

max 


maximum  displacement  of  matrix  in  the  descending  branch 
of  the  uniaxial  tensile  test 


r 
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shear  strength  of  matrix  and  fiber  interface 


effective  crack  over  depth  ratio  in  notched  beam 
aeff 


specimen  ( 
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") 


6 
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FM, 
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*  (Z) 


permanent  deformation 

relative  beam  rotation  of  the  crack  front  in  notched 
beam  specimen 

increments  of  the  compliances  caused  by  the  pressure 
of  force  and  moment,  respectively 
Muskhellshvili ' s  pentential  function 


v 


Poissons'  ratio 
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Uniform  Crack  Closing  Pressure  Applied  on  a  Crack  in  an 
Infinite  Sheet  ^Paris  and  SihJ  . 

Assumed  Uniform  Crack  Closing  Pressure  due  to  Fiber 
Bridging  and  Microcracking  in  a  Cracked  Asbestos-Cement 
Specimen  ^Lenian  and  BunsellJ  . 

Semi-Infinite  Crack  Model  with  an  Assumed  Fiber  Bridging 
Stress  -  o(t)  ^Foote,  Cotterell  and  MaiJ  . 
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Wires  £ Bowling  and  Groves  J  . 
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Hillerborg  J  . 

Fracture  Model  for  Fiber  Reinforced  Concrete. 
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Uniaxial  Tensile  Behavior  of  Fiber  Reinforced  Concrete. 
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Normalized  Post-Cracking  Tensile  Stress-Displacement 
Relationship . 
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Fig.  9  A  Single  Cantilever  Beam  and  the  Applied  Loads. 

Fig.  10  Double  Cantilever  Beam;  Test  Setup. 

Fig.  11  A  Model  for  Fiber  Reinforced  Concrete;  Notched  Beam 

Specimen. 
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Fig.  13  Concept  of  Modified  Strain  Energy  Release  Rate. 
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Fig.  16  Comparison  of  Experimental  and  Theoretical  Load-Line 
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Fig.  17  Comparison  of  Experimental  and  Theoretical  Load-Deforma¬ 

tion  Curves  for  Fiber  Reinforced  Concrete 
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Experimental  and  Theoretical  Compliances  for  Fiber 
Reinforced  Concrete:  (V^  «  IX,  L  ■  0.25",  d  ■  0.006")  _  DCB. 
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Experimental  and  Theoretical  -  Curves  -  DCB  : 
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fig.  7  Uniaxial  Tensile  Behavior  of  Fiber  Reinforced  Concrete 
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Fig.  23  Comparison  of  Theoretical  Permanent  Deformation  for  Different 
Fiber  Volume  Fraction  and  Aspect  Ratio 
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Fig.  24  Comparison  of  Theoretical  Unloading  Compliance  for  Different 
Fiber  Volume  Fraction  and  Aspect  Ratio 


STRAIN  ENERGY  RELEASE,  RATE  .  G0  .  (lb*/ in) 


Fig.  25  Comparison  of  Theoretical  Strain  Energy  Release  Rate  for  Different  Fiber 

Volume  Fraction  and  Aspect  Ratio 
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Fig.  27  Matrix  Process  Zone  and  Fiber  Bridging  Zone  in  Fiber  Reinforced  Concrete  -  DCB  Specimens 


Fig.  28  Matrix  Process  Zone  and  Fiber  Bridging  Zone  in  Fiber  Reinforced  Concrete 

Notched  Beam  Specimens 


